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LONGITUDINAL STABILITY AND CONTROL DERIVATIVES 
OF A JET FIGHTER AIRPLANE EXTRACTED FROM FLIGHT TEST DATA 
BY UTILIZING MAXIMUM LIKELIHOOD ESTIMATION 

By George G. Steinmetz, Russell V. Parrish, 
and Roland L. Bowles 
Langley Research Center 

SUMMARY 

A method of parameter extraction for stability and control derivatives of aircraft 
from flight test data, implementing maximum likelihood estimation, has been developed 
and successfully applied to actual longitudinal flight test data from a modern sophisticated 
jet fighter. The results of this application establish the merits of the estimation tech- 
nique and its computer implementation (allowing full analyst interaction with the program) 
as well as provide data for the validation of a portion of the Langley differential maneu- 
vering simulator (DMS) . The results are presented for all flight test runs in tabular form 
and as time history comparisons between the estimated states and the actual flight test 
data. Comparisons between extracted and manufacturer’s values for five major deriva- 
tives are presented and reveal good agreement for these principal derivatives with the 
exception of the static longitudinal stability derivative C mo ,. This particular derivative 
is extensively investigated by utilizing the interactive capabilities of the computer pro- 
gram. The results of this investigation verify the numbers extracted by maximum likeli- 
hood estimation. 


INTRODUCTION 

A method of parameter extraction for stability and control derivatives of aircraft 
from flight test data has been developed at the Langley Research Center (ref. 1). This 
method utilizes maximum likelihood estimation, an advanced mathematical theory, and 
has been fully implemented on the Langley Research Center real-time digital complex 
(ref. 2). The casting of the digital program in such a framework permits the analyst to 
interact with the program through a control console and cathode ray display. 

Applying the computer program to actual flight test data from a modern sophisticated 
jet fighter airplane (fig. 1) serves a dual purpose: first, to establish the merits of the 
computer program and, second, to provide additional means of validating the Langley 


differential maneuvering simulator (DMS) program (ref. 3). The successful extraction of 
stability derivatives from real-world flight test data goes a step beyond the use of analyti- 
cally generated data toward establishing the creditability and merits of any estimation 
technique. The documentation of the aerodynamic characteristics of the test airplane 
(ref. 4) is extensive and is currently being used as a prime source for the DMS program 
package. The DMS provides dual-cockpit fixed-base simulation in which the two piloted 
vehicles can be operated jointly or individually in either a cooperative mode or an uncoop- 
erative mode. The flight tests themselves and the parameters extracted serve as addi- 
tional means of validating the DMS for this particular airplane. 

The flight test runs utilized in this study are longitudinal responses generated by 
stabilator deflections in the neighborhood of ±5°. The changes in angle of attack and 
pitching velocity are typically ±10° and ±20° per second, respectively. Normal accelera- 
tion responses are between +4g and -lg. The parameters extracted were the standard 
body-axis longitudinal stability and control derivatives. 

SYMBOLS 


Measurements and calculations were made in the U.S. Customary Units. They are 
presented herein in the International System of Units (SI) with the equivalent values in the 
U.S. Customary Units given parenthetically. 
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pitching-moment coefficient 


C m o pitching-moment coefficient at a = 6 = 0 

9C t 


'm„ - 


— per radian 


qc 
,2V ; 

9C m 9C m 

C m "4" C m . damping-in-pitch derivative, — + , per radian 


m q m d 


iw 

9C 


'm 


a 


static longitudinal stability derivative, — — , per radian 


da 


9C 


m • /. _\ 

iw) 


— per radian 


2 


T ,U T 
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longitudinal-force coefficient 
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normal-force coefficient at a = a T , 6 = 6 T 


C y = — — per radian 
^6 96 


mean aerodynamic chord, meters (ft) 

acceleration due to gravity, meters/second 2 (ft/sec 2 ) 

aircraft moment of inertia about the Y-axis, kilogram-meters2 (slug-ft 2 ) 


Mach number 
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mass of fueled airplane, kilograms (slugs) 


pitching angular velocity, radians/second 

pitching acceleration, radians/second 2 

dynamic pressure, i-pV^, newtons/meter2 (lb/ft^) 

Ci 

wing area, meters2 (ft2) 

velocity along horizontal body axis, meters/second (ft/sec) 
acceleration along horizontal body axis, meters/second2 (ft/sec2) 
true airspeed, meters/second (ft/sec) 
velocity along vertical body axis, meters/second (ft/sec) 
acceleration along vertical body axis, meters/second 2 (ft/sec 2 ) 
angle of attack, radians 

rate of change of angle of attack, radians/second 
trim angle of attack, radians 
stabilator deflection angle, radians 
stabilator deflection angle at trim, radians 
pitch angle, radians 

rate of change of pitch angle, radians/second 
mass density of air, kilograms/meterS (slugs/ft^) 

natural frequency of the pitching-moment oscillation for controls locked, 
radians/second 



FLIGHT TESTS 


The flight test data were provided by the U.S. Naval Air Test Center at Patuxent 
River, Maryland. The flight tests were conducted by Navy test pilots as part of an inves- 
tigation with a McDonnell Douglas F-4 airplane. Five different longitudinal response runs 
were made: three during one flight test of the airplane and two during a second flight test. 
The first three runs were made at an altitude of approximately 1524 m (5000 ft) and at 
Mach numbers of about 0.6, 0.7, and 0.8, respectively. The other two runs were made at 
an altitude of approximately 6096 m (20 000 ft) and at Mach numbers of about 0.6 and 0.8, 
respectively. The stability augmentation system (SAS) was deactivated in order to pro- 
vide full response for all the test rims. 

For each of the test runs the procedure was identical. The airplane was trimmed 
by the pilot at the desired altitude and Mach number and held for a short period, and then 
a disturbing input was placed on the stabilators. The pilot task during the maneuver was 
to null any lateral motions. Since thrust was held constant during the maneuver and since 
altitude and Mach number were approximately constant throughout the maneuver, the 
translational modes of the airframe were not sufficiently excited to allow extraction of 
the longitudinal- force coefficient. A time history of the stabilator position for a typical 
test run is shown in figure 2(a). Shown in figure 2(b) are typical time histories of hori- 
zontal velocity u, vertical velocity w, pitching angular velocity q, and normal accelera- 
tion a z . True airspeed was first determined from figure 1 of reference 5 using Mach 
number, pressure altitude, and temperature from flight tests and then resolved through 
angle-of-attack measurements to yield horizontal and vertical velocities. 

Angle of attack was measured by means of a vane on a nose boom. The boom was 
assumed to be rigid. The readings were corrected for pitch by the Navy. Pitching veloc- 
ity was measured by rate gyros located slightly forward and at the foot level of the pilot. 
Normal acceleration was measured by an accelerometer located in the left wheel well. 

The location of the instrumentation is shown in figure 1. However, the center-of- gravity 
offset of the accelerometer was negligible for longitudinal responses. (No documentation 
was available from the Navy as to the accuracy of the instrumentation.) 

The data for each run were sampled at 0.1-sec intervals. Each of these points was 
used by the computer program but every other point is plotted in all time history figures 
to avoid crowding. No smoothing techniques were applied to the data provided. The data 
were supplied to Langley Research Center in tabular form but were transferred to mag- 
netic tape for use with the computer program. 
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COMPUTATIONAL METHOD 


The equations of motion used by the computer program for longitudinal motion are 
as follows: 


u = -g sin 0 - wq + L 


pV 2 S 


2 m 


( C X) Q , T; 6 T + C X a (<* - « T ) + C Xq 


1 PV 2 S 

w = g cos 9 + uq + £ 


2 m 


( C z) a 6 T + C Z a ( a " “T) + % + C Z 6 ( 6 “ 6 T ) 


T 


. I pV 2 Sc 

q = « 


2 L 


( C m) Q!T ,6 T + C m a (« - « T ) + C md |^+ C mq |^+ 0^(6 - 6 T ) 


6 = q 


a z = w - g cos 9 - qu 


a = tan- 12- 
V = ]j u 2 + 

Again, it should be noted that the form of input excitation (pitch doublet yielding the 
short-period motion) makes the extraction of longitudinal-force coefficients ( C X deriva- 
tives) infeasible, since u is negligible. However, the program extracts a set of C x 
derivatives that provide a longitudinal response to match the flight data; thus, a proper 
representation of true airspeed is ensured. 

The method of parameter estimation incorporated into the computer program is 
based on maximum likelihood estimation. This method is inherently iterative and requires 
an initial estimate for each derivative in order to begin the extraction process. Refer- 
ence 1 describes in detail both the method of extraction and the computer program. It 
should be emphasized, however, that the procedure for extracting the derivatives is not 
straightforward. Engineering judgment tempered with estimator statistics, such as the 
variance of individual derivatives and intercorrelation coefficients, plays a significant 
role in the extraction procedure; thus, the capability of analyst program interaction is 
highly desirable. 
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RESULTS 


Figure 3 illustrates the response of the mathematical model to a typical set of 
starting values; the flight data are included for comparison. Usually 8 to 10 iterations 
are required before initial convergence is obtained. From this base, engineering judg- 
ment and statistical information are used to eliminate all possible inconsistencies between 
aerodynamic principles and statistical results. 

The conditions for the five flight tests are listed in table I. In figure 4 are presented 
the model responses generated by the final estimates of the stability derivatives for the 
five flight test runs and the respective flight test data. Again, it should be noted that the 
flight data displayed have not been filtered. 

The total set of extracted stability and control derivatives for each flight test run 
and the standard deviation of each derivative are listed in table Du The intercept deriva- 
tives Cz j0 and C m o are defined at zero angle of attack and zero stabilator deflection. 

From the set of extracted derivatives for the test airplane, the five principal deriv- 
atives Cz , Cz,c, ^m ’ ^m n + Cm ^ , and C m . are plotted as a function of Mach 
number and compared with the manufacturer's data for an earlier version in figures 5 to 
9. (It should be noted that the center-of- gravity variation between the flight test runs is 
less than 0.55% about 31.80% c.) Aerodynamic differences between the two airplanes are 
believed to be small and, hence, two different sets of estimates of the principal derivatives 
exist to predict airplane response. 

The five principal derivatives taken from the manufacturer's data were fixed as 
constants in the program and the minor derivatives ^(Cx) a g ^X a > ^Xq» (C-z)®^ 

C Zq , (Cm)o! T were extracted from the flight data by the maximum likelihood algo- 
rithm. This procedure yielded the model responses presented in figure 10 for the five 
flight test runs. These results clearly indicate an error in both frequency and amplitude 
for normal acceleration a z and an error in amplitude for vertical velocity w. The 
issue now is to uncover the manufacturer's estimates that give rise to the errors. Exam- 
ination of figures 5 to 9 reveals that the largest discrepancy occurs in the estimates of 
the static longitudinal stability derivative C mQ; and, therefore, the error analysis was 
begun with this derivative. 

The principal derivatives taken from the extracted set were used in conjunction with 
the manufacturer's values of C mQ! to generate the responses shown in figure 11; the 
minor derivatives were extracted in the usual manner. These results clearly indicate an 
error in frequency for all the variables. 

In a final attempt to use the manufacturer's value of C mQ! , this derivative was 
fixed at the manufacturer's value and the algorithm was allowed to select the values of 
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all other derivatives. This procedure converged for all flight test runs with the exception 
of test run 4, which diverged. The resulting model responses with the respective flight 
test data are presented in figure 12 for test runs 1, 2, 3, and 5. Again, an error exists 
in both frequency and amplitude for normal acceleration a z and an error exists in ampli- 
tude for vertical velocity w. 

On the basis of these results, it is concluded that the flight test data require a much 
lower value for the static longitudinal stability derivative C mQ! than that presented by 
the manufacturer. As further verification of this conclusion, values of the static longitu- 
dinal stability derivative were calculated by using the following equation: 




“n 2 h 

qSc 


( 1 ) 


These calculations were made from the flight test data after input excitation had essen- 
tially ceased. The results of these calculations are presented in figure 7. The calcula- 
tions yielded values of C m ^ which agreed more closely with the extracted values than 
with the manufacturer's data. 

At this point in the error analysis, at least one discrepancy has been shown to be 
significant. To determine discrepancies among the remaining four principal derivatives, 
the following procedure was used. The four principal derivatives taken from the manu- 
facturer's data were fixed as constants in the program, and the minor derivatives as 
well as CniQ, were extracted in the usual manner. This procedure yielded the model 
responses presented in figure 13 and the variation of C mQ , with Mach number presented 
in figure 14. Although the model responses for this procedure are quite good, the least- 
squares residual for fitting the flight data is inferior to that of the extracted set presented 
in table n. Therefore, it is believed that good agreement exists for all principal deriva- 
tives with the exception of the static longitudinal stability derivative C mQ ,. Also, it 
should be noted that this procedure yielded values for C mQ! that agreed more closely 
with the extracted values than with the manufacturer's data. 


CONCLUDING REMARKS 


It is believed that the agreement between the extracted values and the manufacturer's 
data for four of the five principal derivatives ^ 2 ^, c Z g > Cm q + Cm o!’ C msj> and the 
extensive investigation carried out on the fifth principal derivative Cm^ establish not 
only the merits of maximum likelihood estimation in a real-world environment but also 
the merits of the interactive capabilities of the computer implementation utilized to carry 
out the investigation. The investigation of C mQ , revealed that the flight test data require 
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a much lower value for the static longitudinal stability derivative than that presented by 
the manufacturer. 

The present study also provided some of the data to validate the Langley differential 
maneuvering simulator (DMS) program. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 19, 1972. 
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TABLE I.- FLIGHT TEST CONDITIONS 


Test run 

Altitude 

Mach 

Input 

Excitation 

Center of 

m 

ft 

number 

gravity 

1 

1524 

5 000 

0.6 

Pitch doublet 

Short-period motion 

31.80% c 

2 

1524 

5 000 

.7 

Pitch doublet 

Short-period motion 

31.47% c 

3 

1524 

5 000 

.8 

Pitch doublet 

Short-period motion 

31.27% c 

4 

6096 

20 000 

.6 

Pitch doublet 

Short- period motion 

32.22% c 

5 

6096 

20 000 

.8 

Pitch doublet 

Short- period motion 

31.52% c 


TABLE II.- EXTRACTED STABILITY AND CONTROL DERIVATIVES 


Derivative 

Test run 1 
(c.g. = 31.80% c) 

Test run 2 
(c.g. = 31.47% c) 

Test run 3 
(c.g. = 31.27% c) 

Test run 4 
(c.g. = 32.22% r) 

Test run 5 
(c.g. = 31.52% c) 

Extracted Standard 
value deviation 

Extracted Standard Extracted Standard 
value deviation value deviation 

Extracted 

value 

Standard 

deviation 

Extracted 

value 

Standard 

deviation 

Cm,o 

0.0193 

0.00008 

0.0156 

0.00008 

0.0148 

0.0002 

0.0217 

0.0001 

0.0177 

0.0001 

C m 

lu a 

-.106 

.0017 

-.126 

.0026 

-.167 

.0037 

-.109 

.0009 

-.163 

.0013 

Cmq + CmQ, 

-3.32 

.06 

-3.11 

.08 

-3.30 

.11 

-4.09 

.04 

-3.93 

.05 

c m s 

.576 

.009 

.555 

.015 

.540 

.020 

.618 

.010 

.600 

.010 

c Z,o 

-.080 

.004 

-.086 

.004 

-.106 

.005 

-.025 

.004 

-.049 

.004 

c Z a 

-3.36 

.07 

-3.33 

.10 

-3.64 

.10 

-3.50 

.04 

-3.78 

.06 

C Z6 

.346 

.018 

.334 

.019 

.325 

.022 

.372 

.018 

.361 

.017 

° Z q 




Unidentifiable due to insensitivity 






rad 


Accelerometer 




Time, sec 

(a) Time history of stabilator deflection angle. 
Figure 2.- Typical flight data. 



rad/sec w, ft/sec u, ft/sec 


699 


679 


659 


639 


619 

100 


50 


-50 


— 






4 

l^4t. 







— 


+ ++ + - 

.++++■ 

■' r+ ++- 

+ + + + • 

■+ + + T ‘ 

±111. 

IJ 1 1 

1 LLl 

U LI. 

.1 UJ 

1 1 1 1 

11 IJ 

l“ _ 


+ + 





— 


4 





- 


+ 


, 4 - 

T 

+ 


— 




4 

+ + 






4 

J. 


— 

-f 


4 



- + + +r 

„ 

+ + 


4 . + 




— ■ 




- - 

- 


±111 

JJJJ- 

1111. 

1 1 1 1 

Jill 

J 1 J 1 

1J1L 


l.L 


UJ 




111 




LIU 




'' r 4 , , 

4 ■ 

nn 

t T t 

nn 


213.06 

206.96 

200.86 

194.77 

188.67 

30.48 

15.24 


-15.24 


-30.48 


.50 


4 








I 

.25 

— 


. 






-■ - 

0 

-.25 

- 

4 

4 

4 

.+++ 



■ 4 -r -r , 


_r r r -r 

r 

44 

4 

4. 

4 

4 

.4 

- 

’4 . 
T +r 

,+ + 



+ + + T 

-.50 

nil 

i.lll 

-LIU. 

. 1 Lll 

JILL 

i LJ.LL 

1LLL 

. I.ILJ. 

-LLL) 


11 


I + 

-t + + rf 


1111 


LJ I I 


11 


+ -4-U 


J LI I LUHjlJ.1 1 


■ f+ 


1 111 1 J J II I I I II I 1 II I I 

0 2 4 6 8 10 

Time, sec 

(b) Time histories of horizontal velocity, vertical velocity, 
pitching velocity, and normal acceleration. 


Figure 2.- Concluded. 
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Figure 5.- Variation of C Zq! with Mac ^ number. 


A Extracted from flight; 1524m (5000 ft) 

O Extracted from flight; 6096m (20000 ft) 

...... - Manufacturer's data, c.g.=31%c; 1524 m (5000 ft) 

____ Manufacturer's data, c.g.=31%c; 6096m (20000 ft) 
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Figure 6.- Variation of C Zg with Mach number. 
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Extracted from flight; 1524m (5000 ft) 

Extracted from flight; 6096m (20000 ft) 
Manufacturer's data (extrapolated), c.g.=31 %c; 

Manufacturer's data (extrapolated), c.g.=31%c; 


1524m (5000 ft) 
6096m (20000 ft) 



Figure 8.- Variation of C m q + C m ^ with Mach number. 
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(a) Test run 1. 

Figure 10.- Final model responses generated by using 
manufacturer's derivative estimates. 
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— ■■ Manufacturer's data; 1524 m (5000 ft) 
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□ Algorithm extracted using fixed manufacturer's 
principal derivatives; 1524m (5000 ft) 
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Figure 14.- Variation of C mQ , with Mach number. 


NASA-Langley, 1972 


2 43 


NATIONAL AERONAUTICS AND SPACE ADM I ST RAT I ON 
WASHINGTON, D.C. 20546 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


OFFICIAL BUSINESS 
PENALTY FOR PRIVATE USE J300 


FIRST CLASS MAIL 



007 001 Cl 0 02 720204 S00903DS 

DEPT OF THE AIR FORCE 
A F WEAPONS LAB (AFSC) 

TECH LIBRARY/WLOL/ 

ATTN; E LOU BOWMAN, CHIEF 
KIRTLAND A FB NM 87117 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


"T he aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs., 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




